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ABSTRACT 

We present a timing analysis of the 2009 outburst of the accreting millisecond X- 
ray pulsar Swift J1756. 9-2508, and a re-analysis of the 2007 outburst. The source 
shows a short recurrence time of only ~ 2 years between outbursts. Thanks to the 
approximately 2 year long baseline of data, we can constrain the magnetic field of the 
neutron star to be 0.4 x 10 8 G B H 9 x 10 8 G, which is within the range of typical 
accreting millisecond pulsars. The 2009 timing analysis allows us to put constraints on 
the accretion torque: the spin frequency derivative within the outburst has an upper 
limit of \v\ ^ 3 x 10~ 13 Hzs^ 1 at the 95% confidence level. A study of pulse profiles 
and their evolution during the outburst is analyzed, suggesting a systematic change 
of shape that depends on the outburst phase. 
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1 INTRODUCTION 

The accreting millisecond X-ray pulsar (AMXP) Swift 
J1756. 9-2508 (henceforth referred to as J1756) was discov- 
ered in 2007 by t he Burst Alert Tele scope (BAT) aboard 
the Swift satellite l|Krimm et al.l 12007), Follow up observa- 
tions were performed with the RXTE Proportional Counter 
Array (PCA) which revealed an X-ray pulsar with a spin 
frequency of 182 Hz orbiting in 54.7 min arou nd a very low 
mass companion of minimum mass O.OO67M0 (|Krimm et al.l 
l2007h . 

The X-ray lightcurve of the 2007 outburst showed a very 
fast rise of less than 1 day, followed by an exponential decay 
with an e-folding time of 7.4 ± 0.4 days. An inspection of 
archival RXTE All-Sky Monitor (ASM) data revealed a lack 
of detections in the previous 10 years of observations, which 
suggested a long recurrence time for the outbursts. However, 
due to the low flux of the source and the poor sensitivity 
of the ASM, the recurrence time was only tentative and a 
shorter duty cycle was not firmly excluded. 

In July 2009 a new outburst was disco vered by both 
Swift BAT and RXTE PCA observations (|Patruno et al.l 
12009 j ). Follow-up observations with the Swift X-ray tele- 
scope (XRT) and RXTE PCA were promptly scheduled and 
monitored the outburst of J1756 until it returned into qui- 
escence. 

In this Letter we report a first analysis of the 2009 out- 
burst and a re-analysis of the 2007 data to evaluate the 
long-term evolution of J 1756 on a baseline of approximately 

2 years. 



2 X-RAY OBSERVATIONS AND DATA 
REDUCTION 

We used all the RXTE PCA data (2007 and 2009) and 
the 200 9 Swift-XRT observa tions for J1756 (Table 1). We 
refer to IJahoda et al.l ((2006) for PCA characteristics and 
RXTE absolute timing. We used all available Event and 
GoodXenon data, rebinned to 1/8192 s. For the timing anal- 
ysis, the photons were extracted from the 2.5-16 keV band to 
maximize the signal-to-noise ratio (S/N). We barycentered 
the photons time of arrivals with the FTOOL faxbary by us- 
ing the source posi tions as determined from the Swift-XRT 
2007 observations l|Krimm et al, 2007). We folded approx- 
imately 500 s long data chunks, keeping only those with 
S/N> 3. 

The fractional sinusoidal amplitude of the i-th pulse 
profile is calculated as: 



A, 



N, 



ph,i 



(1) 



where N p h,i and Bi are the total number of photons and the 
background counts (calculated with the FTOOL pcaback- 
est) in the i-th pulse profile, and Ai is the amplitude (in 
photon counts) of the observed pulsations. The error on the 
fractional amplitude Ri is calculated propagating the errors 
on Ai and N p h,i- The error on Bi is negligible with respect 
to the other errors and will not be considered further. The 
X-ray lightcurve is constructed by using the counts in PCA 
Absolute channels 5-37 (~ 2.5 — 16 keV). 

During the 2009 outbursts, the Swift-XST observed the 
faint tail of the outburst with 6 pointed observations (see 
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Table 1. RXTE and Swift-XRT observations analyzed for the 
2009 outburst 



flXTE 


Start 

(MJD) 

55026.00 


End 
(MJD) 
55038.29 


Time 
(ks) 
183 


Observation IDs 

94065-02-01-* 
94065-06-*-* 


SWIFT 


55036.71 


55047.40 


11 


00030952014-00030952017 
00031455001 
00031456001 



Table 1) covering a total of ~ 11 ks. All the data were 
reduced by using the XRT pipeline (v. 0.12.0). Each ob- 
servation lasted between 1 and 2.2 ks. The data were col- 
lected in photon counting (PC) mode and the source and 
background events were extracted in the energy range 0.5- 
10 keV, by using circular regions with radii of 20 arcsec- 
onds. Given the presence of a very bright source in the 
proximity of J1756 (GX 5-1) the first four Swift-XRT point- 
ings were severely c ontaminated by the "stray- light" effect 
|Moretti et al.| [2009) which limited the sensitivity of the ob- 
servations by increasing the background contamination. 

For the spectral analysis we extracted PCA and 
HEXTE spectra from the standard modes following the 
procedures suggested by the RXTE tearrQ, using the lat- 
est available version of the analysis tools by September 1st, 
2009. For the PCA we only used data of PCU 2, which was 
active in all observations. For HEXTE we only used data 
from Cluster B. No energy channels were binned in any case. 
No systematic errors were added to HEXTE data, while we 
used the standard 0.5% for the PCA data. We fitted simul- 
taneously the PCA and HEXTE data in the 3-25 and 25- 
150 keV range, respectively. A single absorbed blackbody or 
a single power law component did not fit well the data, so 
we us ed a standard absorbed wabs (jMorrison fc McCammonl 
1 19831 ) three-components model (disc-blackbody, a power law 
and a Gauss l ine if necessa ry). To fit the data we used 
Xspec vl2.5.0 (lArnaudlll996t ). We tried also to fit compTT 
l|Titarchuklll994 ) and comvPS (|Poutanen fc Svenssonlll996r i 
spectral models. In all cases we obtained a satisfactory fit 
with reduced \ 2 ~ 0.8 - 1.01 for 82 or 83 de grees of freedom 
(depending on the model used) . 



3 RESULTS 

3.1 X-ray lightcurve 

We extracted the 2.5-16 keV X-ray flux from the RXTE- 
PCA and Swtft-XFCT observations. The lightcurve is shown 
in Figure [T] The rising portion of the lightcurve is miss- 
ing because the RXTE monitoring started approximately 
five da ys after the beginning of the outburst (|Patruno et al.l 
l2009af ). The lightcurve shows a slow decay for ~ 5 days be- 
fore entering a fast decay phase (marked with the vertical 
dashed line in the figure) that brings the luminosity down to 
the detection level of RXTE on a timescale of ~ 2 days. The 

1 http:/ /heasarc. gsfc.nasa.gov/docs/xte/xtegof.html 



slow-decay portion of the lightcurve is well fitted by an ex- 
ponential decay law with an e-folding time of 10.2±0.1 days. 
The fast-decay has an e-folding time of 2.7 ± 0.1 days. Both 
these two val u es are longer than the timescales reported by 
iKrimm et all (|2007h for the 2007 outburst (7.6 ± 0.4 and 
0.6 ± 0.3 days for slow and fast decay respectively; see Fig- 
ure [1]). After only 3 days from the last RXTE-PGA detec- 
tion, the Swift-XST was unable to detect the source. The up- 
per limit on the flux level at this time was 100 times smaller 
than that estimated from the last RXTE-PCA observation, 
and indicates that the source was probably returned into 
quiescence. 



3.2 Color-Color Diagram 

We use the 16 s time resolution Standard 2 mode data to 
calculate X-ray colours. Hard and soft colours are defined, 
respectively, as the 9.7-16.0 keV/6.0-9.7 keV and the 3.5-6.0 
keV/2.0-3.5 keV count rate ratios. The energy-channel con- 
version is done by using the pca_e2c_e05v02 table provided 
by the RXTE team. The background was subtracted and 
dead-time corrections were made. In order to correct for the 
gain changes and the differences in effective area between the 
PCUs themselves, we normalized our colours by the corre- 
sponding Crab Nebula co lor values (see lKuulkers et alll994l ; 
van Straaten et"ai1l2003h that are closest in time but in the 
same RXTE gain epoch, that is, with the sam e high volt- 
age setting of the PCUs (jjahoda et al.ll2006h . During the 
2009 observations, the source was always in the so-called 
island/extreme island state, which is typical of persistent 
AMXPs. The source was also in the island/extreme island 
state (not shown) d uring the 2007 outburst, as reported by 
iLinares "et" al. (2008), that identified J1756 as a typical atoll 
source. The colours of the source are consistent among the 
two outbursts. 



3.3 Coherent timing analysis of the 2009 outburst 

Wc fitted each pulse profile with a constant plus two sinu- 
soids: one at the nominal pulse frequency and the other at 
twice the pulse frequency, representing the fundamental and 
the first overtone, respectively. We then selected only pulses 
with a ratio between the amplitude and its statistical error 
larger than 3. Several significant first overtones were found. 
We therefore analyzed separa tely the fundamental a nd the 
first overtone, as described in iPatruno et all (|2009bl ). 

We fitted the time of arrivals of each pulse profile with 
a circular keplerian model plus a linear term representing 
the neutron star spin frequency. We used the pulsar timing 
program TEMP02, which minimizes the timing residuals 
between th e predicted and th e observed arrival times of the 
pulsations (IHobbs et al.ll2006T ). After obtaining a converged 
solution we re-iterated for a few times the whole procedure: 
using the improved solution we re-folded our pulse profiles 
and re-fitted our time of arrivals until convergence. The fit is 
statistically satisfactory for the fundamental (reduced \ 2 = 
1.0 for 251 degrees of freedom) while it is not acceptable for 
the first overtone (reduced \ 2 = 2.8 for 73 dof). The reason 
why the first overtone shows a bad fit relies on the presence 
of X-ray timing noise in the pulse phases. The timing noise 
operates on rather short timescales, producing jumps in the 
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Figure 1. X-ray lightcurve and pulse profile evolution of 2009 (left panel) and 2007 (right panel) outbursts. Left Panel: RXTE-PCA 
[2.5-16 keV] and Swift XRT [0.5-10 keV] lightcurve of the 2009 outburst of J1756. The open circles and the open squares refer to 95% 
confidence level upper limits. The vertical dashed line identifies the beginning of the fast decay. The strong contamination of the nearby 
source GX 5-1 affected both the flux determinations of RXTE-PCA and Swift-XRT (see text for details). The three panels in the center 
show the pulse profile evolution (two cycles are plotted for clarity) obtained when folding a small chunk of lightcurve. Right Panel: 
RXTE-PCA [2.5-16 keV] lightcurve of the 2007 outburst, shown for comparison 




012345678 1 234567 

MJD-55025 [days] MJD-54264 [days] 

Figure 2. Fractional amplitude vs. time for the 2009 (left panel) and 2007 (right panel) outburst of J1756. The fundamental is indicated 
with blue triangles and the first overtone with black filled circles. The fractional amplitude varies between approximately 2 and 20% 
during the two outbursts. The two circled groups of points, that show a substantial increase in pulsed fraction, refer to the beginning of 
the fast decay. 



timing residuals of up to 0.12 cycles on a timescale of a 
few hours. This is not observed in the fundamental, where 
the phase behaviour is smooth throughout the outburst. To 
partially take into account the presence of timing noise in 
the first overtone, we increased the statistical errors on the 
pulse phases of the first overtone by a factor 1.7. The orbital 
and spin parameters of the 2009 outburst, along with upper 
limits on the spin frequency derivative and the eccentricity 
are reported in Table 2. 

We repeated the coherent timing analysis of the RXTE- 
PCA observations of the 2007 outburst of J1756 (see 
iKrimm et al. 1 120071 for a similar timing study of the 2007 
outburst). Also in this outburst several significant first over- 
tones are detected in the pulse profiles. The timing solution 
for the 2007 outburst is reported in Table 3. In this case both 
the fundamental and the first overtone show a low, but sig- 
nificant, amount of X-ray timing noise in the pulse phases. 



We rescaled the pulse phase statistical errors by a factor 1.2 
and 1.4 respectively, to partially take into account this effect 
when estimating the timing parameter uncertainties. 

The pulse fractional amplitudes, as defined in Eq.[T] are 
reported in Figure (note that this is the sinusoidal frac- 
tional amplitude, as defined in Eq. [1] which is a factor \f2 
larger than the fractional rms amplitude). The pulsed frac- 
tion of the fundamental fluctuates between 2 and ~ 20% 
during both outbursts, while the first overtone reaches a 
maximum amplitude of 10-15% and a minimum value of 
2%. The fractional amplitudes increase substantially in cor- 
respondence with the beginning of the fast decay (see circled 
points in Fig. [2] and vertical dashed line in Fig. [T}. Interest- 
ingly, after the beginning of the 2009 fast decay, the frac- 
tional amplitudes return back to values comparable to the 
rest of the outburst. 

We also studied the energy dependence of the pulsed 
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4 DISCUSSION 




10 15 
Energy [keV] 

Figure 3. Fractional amplitude vs. energy for the 2009 outburst 
of J1756. The fundamental is indicated with black triangles and 
the first overtone with black filled circles. The fractional ampli- 
tude of the fundamental increases with energy (apart from a small 
decrease around ~ 6 — 7keV). 



fractions for the 2009 outburst (see Figure [3J). The ampli- 
tude of the fundamental increases with energy, with a slope 
of 0.28±0.02% keV , although it is evident a slight decrease 
in amplitude in the 6-7 keV energy bin. The first overtone in- 
stead remains approximately constant with a slight decrease 
of amplitude above 6 keV. 



4.1 Pulse profiles 

The pulse profile shape of J1756 shows a substantial vari- 
ability during the two outbursts. Pulse shape variability is 
related with the variability of pulse phases and the variation 
of pulse amplitudes with time. The pulse profiles appear to 
have a similar evolution in the two outbursts, with simi- 
lar profile shapes for similar outburst phases (see Fig. [1} . In 
particular, the profiles start with an almost sinusoidal shape 
that evolves in time and turns into a skewed profile during 
the slow decay and becomes double peaked during the fast 
decay stage. Unfortunately, due to low counting statistics, it 
was not possible to further follow the evolution of the pulse 
profiles as the source returns to quiescence. However, it ap- 
25 pears intriguing the s imilarity of this behaviour with that 
of SA X J1808.4-3658 |Hartman et al.ll200Sl . lHartman et"afl 
2009), where the pulse profiles exhibit a similar evolution of 
shapes between outbursts. The presence of a stronger sec- 
ondary bump toward the end of the outburst is also inter- 
esting in light of t he antipodal spot obscu r ation model re- 
cently proposed bv llbragimov fc Poutanenl (|2009l ). The au- 
thors propose a receding inner accretion disc as a possible 
cause of obscuration of the antipodal spot, that slowly re- 
emerges toward the end of the outburst, when the disc is 
emptying. This moves the antipodal hot spot into the line 
of sight of the observer and increases the relative strength 
of the first overtone over the fundamental, thus producing a 
stronger secondary bump in the pulse profile. 



3.4 Emission lines in the X-ray spectrum 

We searched for the presence of emission lines in the broad- 
band spectrum of the longest RXTE observation. In this way 
our statistics are the best and at the same time we are sure 
that the spectrum has not changed significantly in time. 

A line is clearly detected in the 6-7 keV band, in several 
occasions. For example, in Figure [4] we show the ratio of the 
data to the continuum model. In this particular case, the 
flux of the line is (8 ± 2) x 10 -4 photons cm -2 s _1 . We also 
performed a similar analysis on data from the 2007 outburst; 
we found several cases with significant emission lines in the 
same energy range as in 2009. Although we cannot confirm 
the nature of the line in the data, the detection could be 
suggestive of the presence of an Iron line. The low energy 
resolution of our RXTE data does not allow us to constrain 
with high accuracy the energy of the line, the emissivity in- 
dex, the inner radii of the disc and/or the inclination of the 
system. It is also hard to constrain the presence of a Comp- 
ton reflection bump in the X-ray spectrum because of the 
poor statistics given by the very low number of counts at 
high energies. XMM-Newton, Chandra or Suzaku observa- 
tions are needed to verify the nature of the line. 



4-1.1 Timing noise 

The study of the 2009 pulse phases reveals a lack of timing 
noise in the pulse phases of the fundamental, and a short- 
timescale timing noise on the first overtone. This appears 
opposite to previous studies in which the first overtone has 
been considered somehow more stable than the fundamental. 

We repeated the analysis of the 2007 outburst to 
check the long-term spin and orbital evolution of J1756. 
Our orbital s olutio n is consistent with that reported in 
iKrimm et al.l (|2007l ) (see also the erratum lKrimm et al ]|2009l 
for a correction on the value of T a3C ). We found also an offset 
of 0.32 ^tHz from the spin frequency reported in lKrimm et al.l 
(2007). The reason of this discrepancy is more subtle and 
probably has to be identified in the differenc e between the 
timing techniqu es adopted in t h is pap er and in lKrimm et al.l 
l|2007V) . Indeed, IKrimm et "ail (|2007h did not perform an 
analysis of the fundamental and first overtone separately. A s 
shown bv lHartman et al.l (|2008l ) and lPatruno et all (|2009bh . 
significant deviations from the true timing solution arise 
when performing a coherent study in presence of pulse pro- 
file variability. 

The new improved set of ephemeris is not only impor- 
tant for the interpretation of the X-ray data, but is cru- 
cial for understanding the long term evolution of the neu- 
tron star spin and for any radio-search for pulsations in this 
AMXP. 
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Table 2. Timing solution for Swift J1756-2508 (2009 outburst) 



Parameter 


Fundamental 


First overtone 


Spin frequency v (Hz) 


182.06580391(2) Hz 


182.06580393(5) 


Spin frequency derivative i> (Hz/s) 


< 3 X 10~ 13 (95% c.l.) 


< 1.2 x 10 -12 (95%, c.l.) 


Orbital period (s) 


3282.32(3) 


3282.38(6) 


Projected scmimajor axis (lt-ms) 


5.98(2) 


5.96(3) 


Time of passage to the ascending node (MJD) 


55026.03431(3) 


55026.03425(6) 


Eccentricity 


< 0.01 (95% c.l.) 


< 0.01 (95% c.l.) 


X 2 /dof 


250.9/251 


71.0/73 a 



a The value of x 2 refers to the fit where the pulse phase uncertainties have been rescaled by a factor 1.7 



Table 3. Timing solution for Swift J1756-2508 (2007 outburst) 



Parameter 


Fundamental 


First overtone 


Spin frequency v (Hz) 


182.06580393(7) Hz 


182.06580378(24) 


Spin frequency derivative z> (Hz/s) 


< 3 X 10~ 12 (95% c.l.) 


< 1.6 X 10" 11 (95% c.l.) 


Orbital period (s) 


3282.41(15) 


3282.27(26) 


Projected scmimajor axis (lt-ms) 


5.95(4) 


6.04(6) 


Time of passage to the ascending node (MJD) 


54265.28087(10) 


54265.28078(14) 


Eccentricity 


< 0.03 (95% c.l.) 


< 0.03 (95% c.l.) 


X 2 /dof 


59.94/63" 


13.12/14 b 



n The value of x 2 refers to the fit where the pulse phase uncertainties have been rescaled by a factor 1.2. The value of x 2 refers to the fit where the pulse 

phase uncertainties have been rescaled by a factor 1.4. 

4-1.2 Pulse amplitudes 

The pulse profiles of J 1756 have two remarkable features: 
a steep increase of fractional amplitudes with energy and 
a sudden increase of amplitude in correspondence with the 
beginning of the fast decay. 

The energy dependence of the profiles is not new 
among AMXPs, as it has been already obs erved and dis- 
cusse d in three sources: IGR J002 91+5934 dFalanga et al.1 
120051), Aql X-l dCasella et alj|2008l ) and SAX J1748.9-2021 
( Patruno et al.ll2009cj ). A possible reason for this is given by 
photoelectric absorption around the rotating hot spots. A 
very similar slope of the linear relation governing the ampli- 
tude energy dependence w as reported for the i ntermittent 
source SAX J1748.9-2021 (|Patruno et al.1 l2009d V However, 
SAX J1748 and Aql X-l (both interm ittent pulsars) were 
observed to pulsate in the soft state (jCasella et al.1 [20081 . 
lAltamirano et alj |2008), while J 1756 is in the hard state 
(as IGR J00291), which is a more common feature among 
AMXPs. This suggests that the observed increase of pulsed 
amplitudes with energy is not a specific feature of intermit- 
tent sources nor of soft or hard states. 

More difficult to explain is the sudden increase of frac- 
tional amplitude of both harmonics with the beginning of 
the fast decay. Alternate variations of pulse amplitudes af- 
ter the beginning of the fast decay were already observed 
in the pulsed fractio ns of SAX J1808.4-3658 (see Fig. 1 in 
lHartman et~aH 12008ft . These increase might be connected 
with some fundamental readjustment of the accreti on disc 
that corresponds to the onset of the fast decay (see lLasotal 
l200ll for a review of the different outburst stages). 

4.2 Torques during outbursts 

Upper limits on any spin frequency derivative (both spin 
up and spin down) can be made using the coherent timing 
analysis of the 2009 outburst: \u\ < 3 x 10~ 13 Hzs" 1 at the 



Figure 4. Ratio of the data to the continuum model. The flux 



of the line is (8 ± 2) X 10 photons cm 
correspond to ObsID:94065-06-02-04. 



These datasct 



95% confidence level. If standard accretion torque theory 
is correct, the magnetospheric radius (r m ) should decrease 
with the increase of the mass accretion rate M, following 
a power-law r m oc M~ a when r m < r co , with a = 2/7 in 
the simplest case. The term r co is the t he co-rotation radiu s. 
Assuming a source distance of 8 kpc l|Krimm et al.ll2007l ). 
the expected mass transfer rate at the peak of the outburst 
is M ~ 8 x lO~ lo Af0 yr . The expected spin frequency 
derivative is: 



MVGMr„ 



2-kI 



1.6 x lO'^Hzs" 1 



M 



10- 10 M Q yr- 



(2) 
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Therefore the expected maximum spin frequency derivative 
for J1756 has to be of the order of 2 x lCT 13 Hzs" 1 , which 
is within our measured upper limits. 

4.3 Torques during quiescence 

By using the spin frequencies of 2007 and 2009 outbursts we 
can also put some upper limits on the presence of a spin fre- 
quency derivative during quiescence. The difference in fre- 
quencies between the 2007 and 2009 is 0.02/iHz, which is 
within the uncertainties of the frequency determination. By 
propagating the uncertainties on the 2007 and 2009 spin 
frequencies, we can obtain an upper limit for any change 
in spin frequency during quiescence: \/S.v\ ,$ 0.12/iHz (95% 
confidence level) and therefore |z>| < 2 x 10~ 15 Hzs" 1 (95% 
c.L). This value is not particularly constra ining, and is a 
factor 4 larg er than the value measured by iHartman et al.l 
l|20Qg| . 120091 ) for the pulsar SAX J1808.4-3658. 

By using the relativistic forc e-free MHD model s of pul- 
sar magnetospheres described by ISpitkovskvl (|2006T ) , we ex- 
pect a magneto-dipole torque: 

iVdipoic = -V(2Wc) 3 (l + sin 2 a) (3) 

where [i is the magnetic dipole moment and a is the angle 
between the magnetic and rotational poles. This provides an 
upper limit on the magnetic dipole moment and therefore 
on the pulsar magnetic field: 

1/2 



^ < 4.5 x 10 26 (l + sin 2 a 

1/2 



10 45 g cm 2 



(l82 Hz) 



-3/2 



(1 



8 x 10- 15 Hz s- 1 



1/2 



G cm' 



(4) 



By assuming a — 0, we obtain an upper limit for the neutron 
star magnetic field of B < 9 x 10 8 G. If J1756 is located 
close to the Galactic Center at a distance of 8 kpc, a lower 
limit for the magneti c field can be calculated by u sing the 
same argument as in IPsaltis fc Chakrabartvl (|l999h . which 
assumes that the B field is strong enough to channel the 
accreting gas at the maximum mass accretion rate. By using 
the value M « 8 x 1O _1O M0 yr _1 we obtain a magnetic field 
range of: 



0.4 x 10 8 G < B < 9 x 10 8 G 



(5) 



These values are similar to those obtained for other AMXPs 
and confirm the rather narrow range of neutron star mag- 
netic fields in accreting systems. 



4.4 Emission lines and X-ray spectrum 

The presence of an emission line in the 6-7 keV range sug- 
gest the intriguing possibility of an Iron line in the broad- 
band spectrum of J1756. A broad relativistic Iron line has 
been detected only in on e AMXP until now l|Papitto et al.l 
l2009llCackett et al.ll200gh , and its detection in J 1756 can 
help to strongly constrain the inner accretion disc radius 
and the neutron star magnetic field strength. However, care 
has to be used when claiming the detection of Iron lines in 
X-ray binaries, since we do not detect a Compton reflection 
bump in the X-ray spectrum and we do not have a sufficient 



spectral resolution to verify whether the line is composed 
by a superposition of narrow lines or whether it is a single 
genuinely relativistically broadened line. The reason why we 
do not detect a Compton bump in our data can be due to a 
lack of statistics at high energies, and a deeper observation 
is required to verify the Iron line hypothesis. 



5 CONCLUSIONS 

Swift J1756. 9-2508 has shown a short recurrence time of only 
two years and it is the optimal candidate to study the long- 
term evolution of neutron star spin and orbital parameters. 

The pulse profiles evolve in time with strong similari- 
ties between the two outbursts and show a steep increase 
of pulsed fraction with energy. With the increasing number 
of AMXPs exhibiting an hard energy dependence of pulse 
amplitudes, it becomes more likely to discover new AMXPs 
among the large family of the non pulsating low mass X-ray 
binaries. Indeed, all the 12 AMXPs discovered until now 
have been observed with X-ray telescopes which are sensi- 
tive to energies below 20 keV, or that are not designed to 
perform efficient coherent X-ray timing observations. With 
the forthcoming X-ray telescopes like ASTROSAT, which 
will have a very high effective area up to ~ 80 keV, the 
detection of AMXPs with rising pulsed fractions at hard 
energies could be much easier and efficient. 

The pulse phases show a low, but significant, timing 
noise content and special techniques have to be applied to 
determine the neutron star spin evolution. Thanks to the 
measurement of timing parameters in two outbursts, we can 
now set the first constraints on the neutron star magnetic 
field (0.4 x 10 8 G < B < 9 x 10 8 G) that makes this neutron 
star a typical member of the AMXP family. 
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